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EXPONENTIAL TREND TO EQUILIBRIUM FOR THE INELASTIC BOLTZMANN 
EQUATION DRIVEN BY A PARTICLE BATH 


JOSE A. CANIZO & BERTRAND LODS 


Abstract. We consider the spatially homogeneous Boltzmann equation for inelastic hard 
spheres (with constant restitution coefficient a G (0, 1)) under the thermalization induced 
by a host medium with a fixed Maxwellian distribution. We prove that the solution to the 
associated initial-value problem converges exponentially fast towards the unique equilib¬ 
rium solution. The proof combines a careful spectral analysis of the linearised semigroup 
as well as entropy estimates. The trend towards equilibrium holds in the weakly inelastic 
regime in which a is close to 1, and the rate of convergence is explicit and depends solely 
on the spectral gap of the elastic linearised collision operator. 
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1. Introduction 

We pursue our investigation initiated in [3, 4] of the qualitative properties of inelastic 
hard spheres suspended in a thermal medium. In a more precise way, we investigate here 
the large time behavior of the one-particle distribution function f{v, t), v G t > 0 
solution to the following spatially homogeneous Boltzmann equation: 

dtf = QMJ) + Uf), ( 1 - 1 ) 

where Qa{f, /) is the inelastic quadratic Boltzmann collision operator, while L(/) models 
the forcing term. The parameter a is the restitution coefficient, expressing the degree of 
inelasticity of binary collisions between grains: 0 < a ^ 1, and the purely elastic case is 
recovered when a = 1. 

1.1. Setting of the problem. As is well documented, dilute granular flows can be de¬ 
scribed by kinetic models associated to suitable modifications of the Boltzmann operator 
for which hard-sphere collisions are assumed to be inelastic [6] : each encounter dissipates 
a fraction of the kinetic energy. In absence of energy supply the system cools down and 
the corresponding dissipative Boltzmann equation admits only trivial equilibria. This is no 
longer the case if the spheres are forced to interact with an external thermostat, in which 
case the energy supply may lead to a non trivial steady state. Different kinds of forcing 
term have been considered in the literature [10, 21, 15, 16, 17, 11], and the qualitative 
properties of the corresponding steady states have been investigated. In particular, the 
following questions have been addressed regarding steady solutions of kinetic equations 
of the type (1.1): (1) Existence [15, 11], (2) Uniqueness in some weakly inelastic regime 
corresponding to a close to 1 [16, 17] and (3) stability, i.e. convergence of the solution 
to the associated Boltzmann equation towards the steady state [16, 17]. In particular, for 
hard spheres subject to diffuse forcing, the large-time behaviour of the solution to the BE 
has been completely characterised in [17] whereas, for anti-drift forcing (closely related 
to self-similar solutions to the freely evolving Boltzmann equation), asymptotic behaviour 
has been considered in [16]. 

In this paper we are concerned with the asymptotic behaviour of a physical model in 
which the system of inelastic hard spheres is immersed in a thermal bath of particles at 
equilibrium, already investigated in [3, 4] . In this model the forcing term is given by a 
linear scattering operator describing elastic collisions with the background medium: 

Uf) = Qi{f,M) 

where M. stands for the distribution function of the host fluid, supposed to be a Maxwellian 
with unit mass, bulk velocity uo G and temperature 0o > 0: 

The precise definitions of the collision operators Qa{f,f) and L(/) are given in Subsec¬ 
tion 2.1. We refer to [15, 16, 17, 11] for a mathematical discussion of various models and 
their physical motivation. We restrict our attention to interacting hard spheres and refer 
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to [7] for exhaustive references on the pseudo-Maxwell approximation. A salient feature 
of both collision operators Qaif, /) and L is that their only collision invariant is mass, i.e. 



In contrast with the elastic Boltzmann operator, neither the momentum nor the energy are 
conserved by or L. 

The existence of smooth stationary solutions for the inelastic Boltzmann equation under 
the above thermalization has already been proved in [3], for any choice of the restitution 
coefficient a. This has been achieved by controlling the Lt’-norms, the moments and the 
regularity of the solutions for the Cauchy problem, together with a dynamical argument 
based on the Tychonoff fixed-point theorem. 

Uniqueness of the steady state is proven in some previous contribution [4] for a smaller 
range of parameters a. Namely, the main results of both [3] and [4] can be summarized 
as follows: 

Theorem 1.1 (Existence and Uniqueness of the steady state). For any g ^ 0 and a G (0,1], 
there exists a steady solution Fa G L 2 (M^), Fa{v) ^ 0 to the problem 


Qa{Fa^ Fa) + — 0 


(1.3) 



Moreover, there exists ao G (0,1] such that such a solution is unique for a G {ao, 1]. The 
(unique) steady state Fa for a G {ao, 1] is radially symmetric and belongs to C°°(]R^). 

We assume in the sequel that a G (oq, 1] and we are interested in the large-time be¬ 
haviour of the solution f{t,v) to (1.1) (whose existence and uniqueness are guaranteed 
by [3]). In particular, if one assumes 



(1.4) 


then one expects the solution f{t, v) to converge towards the unique steady solution with 
unit mass given in Theorem 1.1. Our goal in the present paper is to provide sufficient 
conditions on the initial datum /o ensuring that this convergence holds true, with an ex¬ 
ponential rate that we make explicit. As in [16, 17], this exponential trend to equilibrium 
is proven to hold in the weakly inelastic regime, i.e. for a range of parameters a G {af 1] 
for a certain explicit > oq. 

1.2. Main result and strategy of proof. Our goal is to prove a quantitative version of 
the return to equilibrium for the solution to (1.1). Our main result combines local sta¬ 
bility estimates in a certain weighted -space with suitable entropy estimates. A crucial 
point in our approach is that it strongly relies on the understanding of the elastic problem 
corresponding to a = 1. In the elastic case, as is well known, Fi = Ad is exactly the host- 
medium Maxwellian appearing in L. The spectral properties of the linearised operator 
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around A4 have been studied in [4]. Namely, in the weighted space 

A'= exp(a|r;|) dti), a > 0, 

the elastic linearised operator given by 

Jfih = Qi(h,A4) + Qi(j\4, h) + L/i, h e (M^; (1 + lul) exp(a|t>|) du) 

admits a positive spectral gap i/ > 0 which can be explicitly estimated. We shall consider 
solutions to (1.1) associated to a nonnegative initial datum /o G satisfying (1.4) and 

-f^(/o|-A4) = / foiv) log ( ) d?; < oo. (1.5) 

JiRS 

Our main result can be stated as follows: 

Theorem 1.2. For any 0 < < n (with v equal to the size of the spectral gap of Jtfi) there 

exists some explicit G (0,1) such that, for any a G (a^, 1] and any nonnegative initial 
datum with /o G -T satisfying (1.4) and (1.5) the solution f{t,v) to (1.1) satisfies 

\\f{t)-Fa\\x^Kexp{-v^t) Vt ^ 0 
for some positive constant K depending on e, a and H{fo\M.). 

Notice that the rate of convergence is explicitly computable in terms of the spectral gap 
of the elastic linearised operator in X. 

The proof of the above is based upon two main ingredients: 

(1) A local stability estimate in which exponential convergence is established for small 
perturbations of the equilibrium state, i.e. whenever the initial datum /o is close 
enough to Fa ■ 

(2) Suitable entropy estimates as a tool to pass from local to global stability. We take 
advantage of the fact that the scattering operator L is dominant in the weakly 
inelastic regime. 

To tackle the above first point (1), we have to perform a fine study of the spectral pro¬ 
perties of both the linearised operator Jtfa around the steady solution Fa and its associated 
evolution semigroup. More precisely, introduce 

^ah = Qa{h,Fa) + Qa{Fa,h) + Lh, h G L^(M^ ; (1 + |r'|)exp(a|tj|) dv), a G (ao, !]• 

We deduce the spectral properties of ^a in X from those of the elastic operator by a 
perturbation argument valid for a close enough to 1. Notice that the elastic limit a —^ 1 
is actually well behaved since the operator gap (in the sense of [13]; see Appendix A) 
between ££a and is going to 0 as a —> 1. This allows us to apply results from the 
perturbation theory of unbounded operators [13] . This strongly contrasts with the analysis 
of [16, 17] which, though perturbative, was ill-behaved in the elastic limit. 

As is well known, the spectral properties of the Co-semigroup {Sa{t))t^o generated by 
J^’a cannot be directly deduced from those of because of the lack of spectral map¬ 
ping theorem in infinite dimensional Banach spaces. In particular, one cannot directly 
derive from the existence of a spectral gap for ^a the decay of the associated semigroup. 
However, following an operator splitting strategy introduced in [18], we can localise the 
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essential spectrum of {Sa{t))t^o through a weak compactness argument and deduce from 
that the local stability theorem (see Theorem 3.9 and Theorem 4.2). We give a direct and 
elementary proof that does not rely on the recent results of [12, 18] . 

In order to address the above point (2) entropy estimates play a crucial role. Our 
method is based upon the following entropy-entropy production estimate recently ob¬ 
tained in [5]. Introducing the entropy production associated to L, 

the result reads as follows: 

Theorem 1.3. There exists A > 0 such that 

D(/) > AiJ(/|A,) := A/W log do 0 

holds for any probability distribution f G L^(R^, dn). 

This result has important consequences on the asymptotic behaviour of the solution to 
(1.1) in the elastic case a = 1. In this case the scattering operator L becomes dominant 
and forces the solution f{t,v) to (1.1) to converge exponentially fast towards the unique 
equilibrium state of L, which is the Maxwellian M. (see [5] for details). Roughly speaking, 
in the elastic limit a —)• 1, we expect the persistence of this behaviour and we expect the 
scattering operator L to drive the system in some neighbourhood of M. Since TA ~ Ad 
for a close to 1, the dynamics is forced to take the solution close to as t —oo. To be 
more precise, using the above Theorem 1.3, one can estimatethe evolution of the relative 
entropy along the solutions to (1.1) (see Proposition 5.1) to get 

H{f{t)\M) ^ exp{-Xt)H{fo\M) + K{1 - a) Vf ^ 0 ; a G (oq, 1] 

for some positive constant K > 0 independent of a. The above estimate ensures that, for 
large time and a ~ 1, the solution to (1.1) will become close enough to M and, hence to 
Fa which, combined with the local stability theorem, yields our main result. It is worth 
mentioning here that, while the analysis in [4] dealt with a (possibly) inelastic scattering 
operator, we restrict ourselves here to elastic interactions between the hard spheres and 
the host medium due to the inavailability of Theorem 1.3 in the inelastic case. Notice 
however that all the spectral results of the paper, as well as the local stability theorem 4.2, 
hold true without modification substituting L by the inelastic scattering operator 

Le/=Qe(/,Af), 

associated to a general constant restitution coefficient e €(0,1]. 

1.3. Plan of the paper. The organisation of the paper is as follows. In the next section 
we define the collision operators Qa and L, and recall from [4] the main properties of 
the solution to (1.1) and the steady state Fa- Section 3 is devoted to a study of the 
spectral properties of and {Sait))t^o, used later in Section 4 to derive the local stability 
Theorem 4.2. In Section 5 we exploit the entropy estimates to establish our main global 
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Stability result. The proof that, for any a E (oq, 1], generates a Co-semigroup in X is 
postponed to Appendix A and uses a weak compactness argument. 

1.4. Notation. Given two Banach spaces X and Y, we denote by ^{X,Y) the set of 
linear bounded operators from X to y and by || • \\^(x,y) associated operator norm. If 
X = Y,we simply denote ^{X) := ^^{X,X). We denote then by ^{X) the set of closed, 
densely defined linear operators on X and hy JY{X) the set of all compact operators in 
X. For A € ^{X), we write ^(A) c X for the domain of A, Yy'(A) for the null space 
of A and Range(A) c X for the range of A. The spectrum of A is then denoted by S(A) 
and the resolvent set is p(A). For A E p(A), R{\,A) denotes the resolvent of A. We also 
define the discrete spectrum 6d(^) as the set of eigenvalues of A with finite algebraic 
multiplicity (see [13, 9] for more details). We denote by s(A) the spectral bound of A, i.e. 

s(A) = supjRe A; A E ©(A)}. 

There are several definitions of the essential spectrum of A in the literature which are un¬ 
fortunately not equivalent. In the present paper we adopt the notion of Schechter essential 
spectrum, denoted by ©ess (A) and defined by 

®ess (•4)= n ©(A +a:). 

KejP{X) 

For a bounded operator T & SS{X) we can also define the essential radius of T as 
ress(r) =inf{r >0;©(r)n{AEC; |A| > r} C ©d(T) } 

= sup{|/i| ; p E &ess{T)} . 

Notice that the first identity is peculiar to Schechter essential spectrum whereas the second 
one is valid for any of the various notions of essential spectrum (see [8, Corollary 4.11, p. 
44]. If iU{t))t^o is a Co-semigroup in X with generator A, we denote by (jJo{U) its growth 
bound and by uiessiU) its essential type, defined by 

exp(ta;o((/)) = sup{|/i| ; p E ©(C/(f))} , 

exp (t (Jess (C)) = ress(C(t)) = sup{|/i| ; p E ©ess(C(f))} fort ^ 0. 

2. Preliminary results 

2.1. The kinetic model. Given a constant restitution coefficient a E (0, 1), one defines 
the bilinear Boltzman operator Qa for inelastic interactions and hard spheres by its action 
on test functions V’(r'): 

[ Qa{f,g){v)^p{v)dv= [ [ [ fiv)g{w)\v-w\(ijj{v')-'ilj{v))dvdwda (2.1) 

with v' = v + {\v — w\a — V+ w). In particular, for any test function one has 

the following weak form of the quadratic collision operator: 

/ QMJ){v)'ip{v)dv = l-[ f f{v) f{w)\v-w\Aa['il^]{v,w)dwdv, 

JR3 ^ Jr3 Jr3 


( 2 . 2 ) 
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where 


AaW\{vjw) = ^ + - i>{v) -'il){w))6.a 

47r 


/S2 (2.3) 

= -4+ [v^] {v,w)- A~ i^p] {v, w) 

(where we have used the symmetry of the integral under interchange of v and w) and the 
post-collisional velocities {v', w') are given by 

1 “I" o 


v' = V + {\qW - q), 


w = w — 


{\qW-q), q = v-w. 


(2.4) 


4 4 

In the same way, one defines the linear scattering operator L by its action on test functions; 

1 f 


Mf){v)'<P{v)dv = 


f{v)M.{w) I?; — w\J[ip]{v, w) drc dt;, 


where 


JbP]{v,w) = ^ [ (ipiv*)-tpiv)) da = J^[tp]{v,w) - J ['Ip]{v,w). 
4vr y§2 


(2.5) 


( 2 . 6 ) 


with post-collisional velocities {v*, w*) 

V* = V+ ]^{\q\a - q), w* = w-]^{\q\a - q), q = v - w. (2.7) 

For simplicity, we shall assume in the paper that the particles governed by / and those 
with distribution function M. share the same mass. Notice that 

Uf) = Qi{f,M) 

and we shall adopt the convention that post (or pre-) collisional velocities associated to 
the coefficient a are denoted with a prime, while those associated to elastic collision 
are denoted with a *. We are interested in the large time behaviour of solutions to the 
following Boltzmann equation: 

= fiO,v) = fo{v), t>0,vGR^. (2.8) 

Notice that 

QMJ) = Qt if, f)- Qaif,f) = Qtif,f) - mf) 

where 

^(/)(^) = (/ * I • l)(^) = / fiw)\v - w\ dw. 

7k3 

Notice that S(/) does not depend on the restitution coefficient a G (0,1]. In the same 
way, 

Mf){v) = L+if){v) - L-if){v) = L+if){v) - ^{M)iv)fiv) 

Existence and uniqueness of solutions to (2.8) have been established in [3]. In particular, 
if /o is a nonnegative initial datum with 

fo{v)\v\^ dv < QO 


j 


and 


fo{v) dv = 1 


(2.9) 
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then there exists a unique nonnegative solution {f{t,v))t^o to (2.8) which additionally 
satisfies 

sup / /(t, f )|r;|^ du < oo, and / f{t,v)dv = l Vt^O. 

JR3 Jr3 

More generally, uniform propagation of moments holds: namely, for any p ^ 2 one has 


[ fo{v)\v\^ dr; < oo 
JR3 



/(t, v)\v\^ dt) < oo. 


( 2 . 10 ) 


See [3, Proposition 4.2 & Theorem 4.8] for more details. Owing to the above mass con¬ 
servation property we shall restrict ourselves in the sequel to nonnegative initial data 
satisfying (2.9). Due to the influence of the scattering operator L there is no additional 
conservation law besides mass conservation. In fact, it appears impossible to express the 
evolution of the momentum 

= / f{t,v)vdv£R^ 
iR3 


and the energy 


in a closed form. 


E{t)= j f{t,v)\v\^dv 
JR3 


2.2. A posteriori estimates. We collect here several results obtained in our previous con¬ 
tribution [4] regarding the properties of solutions to (2.8) as well as those of the steady 
solution Fa to (1.3). We begin with high-energy tails for the solution to (1.1) and Fa- 

Theorem 2.1. Let fo be a nonnegative velocity distribution with fo{v)dv = 1. Assume 
that fo has an exponential tail of order s G (0,2], i.e. there exists tq > 0 and s G (0,2] such 
that 

/ fo{v) exp (rolt'l'') du < oo. 

JR3 

Then there exist 0 < r ^ ro and C > 0 (independent of a ^ (0, l]j such that the solution 
{f{t,v))t^o to the Boltzmann equation (2.8) satisfies 


sup / /(t, v) exp (r|r;|®) dt; ^ C < oo. (2.11) 

t^o Jr3 


In particular, there exist constants A > 0 and M > 0 such that for all a G (0,1] and all 
solutions Fa to (1.3) one has 



exp (A|?;|^) df ^ M. 


Notice that the above integral tail estimate for Fa can actually be strengthened to get 
the following pointwise Maxwellian bounds: 


Proposition 2.2 ([4, Theorems 4.4 & 4.7]). There exist two Maxwellian distributions M 
and M. (independent of a) such that 

M_{v) ^ Fa{v) ^ Al(f) Vf G 


Vo G (oq, 1). 
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2.3. Convergence of Fa to Let us now introduce 

X = {m~^) = , m~^{v) dv), y = L\{m~^) = L^{R^,{v)m~^{v)dv) (2.12) 

where 

(t>) = (1 + lup) 2 and m(t>) = exp (—alul), a > 0, ?; € M^. 

According to the above Theorem 2.1, Fa ^ X for any a G (0,1]. We recall from [1, 
Proposition 11] that Qa is well defined on y\ 

Proposition 2.3. There exists C > 0 such that, for any a G (0,1) 

\\Qaih,g)\\x + \\Qaig,h)\\x ^ C\\h\\y ll^lly V/l, 5 G 3^. 

Moreover, one has the following: 

Proposition 2.4 ([16, Proposition 3.2]). For any a,a' G (0,1) and any f G 
g G L\{m~^), it holds 

\\Qt{f,9) - Qa'if^9)\\x ^p{a-a')\\f\\wf\m-^) l^b 

and 

WQaidJ) - Qa'{9,f)\\x ^P(a-a')ll/llwl’l(m-l) 

where p{r) is an explicit polynomial function with lim,,^o+ P(^) = 0. 

In the elastic limit a —)• 1, one has the following: 

Theorem 2.5 ([4, Theorem 5.5]). There exists an explicit function gi{a) such that limc^i gi{a) 
0 and such that 

Iba --Adb ^ VaG(ao,l]- 

2.4. Spectral properties of the linearised operator for a = 1. Define the elastic lin¬ 
earised operator SF\ : ^(Afi) c A —)■ A by 

= Qi(M,h) + Qi(h,M) + Lh, V/i G ^(ifi) = 3^. 

(We recall X and y were defined in (2.12).) We introduce also 

x = {feX; [ fdv = 0}, y = {fey-, [ fdv = o}. 

3r3 Jr3 

One has the following structure of the spectrum of Ffi: 

Theorem 2.6 ([4, Theorem 5.3]). The null space o/Afi in X is given by 

c/K(.ifi) = span(Al). 

Moreover, SF\ admits a positive spectral gap u > 0. In particular, cA'(.ifi) nX = {0} and 
is invertible from y to X. 
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Let us spend a few words on the strategy used to prove the above result since we will 
use several of the tools involved to study the spectral properties of the linearised semi¬ 
group in the next section. The proof of the above result is related to a general strategy 
introduced in [12] which consists in deducing the spectral properties in from the much 
easier spectral analysis in L?. The existence of a spectral gap for the linearised collision 
operator in ^ is relatively easy to obtain through a suitable Poincare-like 

inequality and the task is to prove that the linearised collision operator in X can be de¬ 
duced from the one in the Hilbert setting. This is done thanks to a suitable splitting of the 
linearised operator as 

^i=Ai+ Bi 

where 

{i) Ai : X ^T-L is bounded; 

(ii) the operator Bi : ^{Bi) —^ X (with ^{Bi) = y) is /3-dissipative for some positive 
/3 > 0, i.e. 


/ 


signf{v)Bif{v)m (t>) dt> ^-/3||/||y V/G 3^. 


(2.13) 


Under these conditions [12] asserts that the spectrum of .ifi in X will be the same of that 
in 

To be more precise, the splitting is as follows. Let us introduce the linearised Boltzmann 
operator 

Tif = Qi{MJ) + Qi{f,M) 

so that .ifi = Tl + L. Clearly, 


Tif = T,+f- 


'i^f-(^i{v)f{v)-M{v) [ f{w)\v-w\dw, 

Jm.3 


while 


Lf{v) = L+{f)-'S{v)f{v) 


and the splitting consists in setting, for some R > 0 large enough so that (2.13) holds, 

Ai= A\+ Af 

with 

•^i/=+ L+(xnfl/), A‘if{v) = -M{v) [ f{w)\v - w\dw 

Jm.3 

where is the open ball in with radius R > 0 and center 0. Then, simply sets 
Bi = y’l — Ai- Notice that, in the above inequality (2.13), the constant /3 > 0 can be 
chosen as 

/) = ^ + Oj_ + E 

for some arbitrarily small e > 0 where 


E = inf XA > 0 


1 + P 

Notice fi does not depend on R. 


and 


= i„f JTW > 0. 

- 1 + f 
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3. Spectral analysis of the linearised operator and its associated semigroup 

We recall that for a G {ao, 1], Fa denotes the unique steady state with unit mass, 
solution to (1.3). In order to study the stability of Fa for a close to 1 we will first prove 
that the following operator has a spectral gap in X ; 

^a{h):=Qa{h,Fa) + Qa{Fa,h) + Uh) {h^y). (3.1) 

Notice that thanks to Proposition 2.3 the above expression is well defined and belongs to 
X whenever /i G 3^. It is fundamental here that the domain of ^a in X does not depend on 
a, i.e. 

^(Jfa) = ^yfi) = y VaG(ao,l]. 

This is in major contrast with the situations investigated in [16, 17] where the forcing 
term was a differential operator for which the domain of the associated linearised operator 
involved Sobolev norms. 

Since Jfa is the linearisation of the nonlinear operator Qaif, f) + L(/) near its steady 
state Fa, a study of its spectrum will allow us to perform a perturbative study of the 
evolution equation (2.8). We deduce from the results of Section 2.3 the following technical 
result stating that J^’a is close to Ffi for a close to 1. It will play a crucial role in our 
analysis: 

Proposition 3.1. There exists an explicit function zu: (ao, 1] —M+ such that lim„^i+ Tu{a) = 
0 and 

\\^a{h) - ^l{h)\\x ^ w{a) ||/l||y V/l G 3^. 

Proof A fundamental observation is that the domain of Ffa is actually independent of a, 
i.e. 3>{^a) = y for any a G (oq, 1]. Let /i G T’ be fixed. We have 

\\^a{h) - ^l{h)\\x = \\Qa{Fa,h) + Q«(/l,F„) - Ql{M,h) - Ql{h,M)\\x 

^ \\Qa{Fa — -M., h) + Qa{h, Fa — A1)||;y 

+ \\Qa{M,h) - Qi{M,h)\\x + \\Qa{h,M) - Qi{h,M)\\x. 

Thus, using Proposition 2.3 and Theorem 2.5, one sees that there exists some constant 
C > 0 and some explicit function ryo(a) with lim^^i po{a) = 0 such that 

WQaiFa -M,h) + Qaih,Fa - M)\\x ^ C\\Fa - M\\y\\h\\y ^ Vo{a)\\h\\y. 

In the same way, according to Proposition 2.4, 

\\Qa{M,h) - Qi{M,h)\\x + \\Qa{h,M) - Qi{h,M)\\x ^ Vi{a)\\h\\y 


some explicit function pi{a) with lima^i r/i(a) = 0. These two estimates give the result 

with ro(-) = %(•)+ n 
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3.1. Spectral gap of the linearised operator. We investigate here the spectral properties 
of in X. We begin by studying the kernel of 

Proposition 3.2. There exists some explicit ai E {ao, 1] such that, given a E (ai, 1], 0 is a 
simple and isolated eigenvalue of atid there exists Ga ^ y with unit mass and such that 

^{yfa) = Span(GQ,) Va E (ai, 1]. 

We denote then by Pq the spectral projection associated to the zero eigenvalue ofj^’a. Then, for 
any f G X, one has Pq,/ = pf Ga with pf := fj^g f{v) di;. In particular, Range(I — Pq) = X 
for any a E (ai, 1]. 

Proof For any a E (oq) !]> set for simplicity Tq, = .ifi — with domain ^{Ta) = y. From 
Proposition 3.1, 

\\T^h\Ui^w{a)\\h\\y \fhGl^{^i)=y. 

Since || • ||y is equivalent to the graph norm of there exists c > 0 such that 

ll^b^c(||h||;, + ||.5^lh||;t) VhE3^ 
from which the above inequality reads 

\\TMx ^ a\\h\\x + b\\^ih\\x Vh E 

with b = cw{a). Since limQ,^i+ w{a) = 0, this makes Tq, a .ifi-bounded operator with 
relative bound b < 1 for any a E (ao,l] for some explicit E (ao,l]- In particular, 
according to [13, Theorem 2.14, p. 203] (cf. Theorem A.2), the gap 6{yfa,y^i) between 
.ifo = yfi + Ta and (as defined in [13, IV.2.4, p. 201]; see Appendix A) is less than 
. Now, recall that the spectrum of .ifi splits as 

e(.ifi) = { 0 } u 6'(.i^i) 

where ©'(.ifi) c {z E C; Rez ^ Denoting by Pi the spectral projection asso¬ 

ciated to the 0 eigenvalue, one gets that X = X (B X with X = Range(Pi) and 
X' = Range (I — Pi) with moreover 6 (.ifi|^//) = {0} and ©(.ifi|^') = 6 '(ifi). In particu¬ 
lar, the two above parts of ©(ifi) are separated by the closed curve 7 ^ = { 2 : E C; \z\ = r}, 
for any r E (0, u). Then, according to [13, Theorem 3.16, p. 212 & IV.3.5] (see Theorem 
A.3), there exists 5 > 0 such that the same separation of the spectrum and decomposition 
of X hold for any operator 5 E "^(W) for which the gap 5(5,ifi) < 6. Choosing now 
ag C («()) 1] such that 5(ifQ,ifi) < 5 as soon as a E (a'g', 1], one gets therefore that the 
spectrum of ©(ifa) can be separated by 7 ^, i.e. it splits as 

e{^a) = e"(ifQ) U ©'(i^Q) Va E (a'o', 1 ] 

where ©"(ifa) C {2 E C; |2;| < r} while ©'(ifa) C {z E C; |2:| > r}. Moreover, the space 
X splits as X = X" © X' with ©(ifaW//) = ©"(ifa) and Moreover, 

still using Theorem A.3, dim(A’Q) = dim(A’Q) = 1 . This shows that actually 

©"(ifa) = {/Xq} 
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where is a simple eigenvalue of ££a with \p,a\ < r for any a G (oq, 1]. Let us show that 
actually pa = Q (at least for sufficiently large a) Define Pq as the spectral projection 
operator associated to pa, i-e. 

Pa = ^ jf VaG(ai,l]. 

According to A.3, one also has liniQ^i ||Pq — Pi = 0 with an explicit rate, from which 
there exists some explicit ai G (og, 1] such that 

llPaZ-Pi/k < 1 forallaG /G (3.2) 


Let us prove that /Xq, = 0 for any a G (oi, 1]. Let us argue by contradiction and assume 
there exists a G (oi, 1] for which pa / 0. Let (pa be some normalized eigenfunction of ^a 
associated to pa, i.e. pa ^ y\ {0} satisfies ^faPa = fJ-a Pa- Integrating over we get that 

[ pa{v) dv = 0. 

Jk3 

For any / G A", there exists /3 = I3{a,f) such that Pq,/ = (3pa while Pi/ = gjAi. In 
particular, one sees that 

[ Fafdv = 0 while [ Fif dv = Qf V/ G 

JK3 Jr3 

This clearly contradicts (3.2). Therefore, for any a G (ai,l], Pa = 0 and the above 
reasoning shows that any associated eigenfunction pa is such that 

[ pa{v) dv^O. 

Let then Ga he the unique eigenfunction of ££a associated to the 0 eigenvalue with 
/r 3 G'q(i’) dv = 1. From [13, Eq. (6.34), p. 180], one has Range(I — Pa) C Range(.ifa) for 
any a G ( 02 ,1]. Since ^af du = 0 for any / G ^(Afa) we get that Range(I — Pa) C Af 
for any 02 < a ^ 1- Thus, given f G X, since / = Pa/ + (I — Pa)/, it holds 


■= f /du = / ^afdv. 

JR3 Jr3 

Since moreover Pa/ = PfGa for some /3j G M and Ga is normalised, we get that Pa/ = 
QfFa- In particular, if / G A" then Pa/ = 0 and / G Range(I — Pa) which achieves the 
proof of the result. □ 


From now on, u > 0 will denote the size of the spectral gap of Afi (see Theorem 2.6). 
Our main result in this subsection is the following: 


^Notice that this cannot be deduced directly from the fact that r > 0 can be chosen arbitrarily small 
since the range of parameters (qq, 1] for which the above splitting holds actually depends on r through the 
parameter S in A.3. 
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Theorem 3.3. Take 0 < l-* < There is ao < < 1, with a* depending on v^, such that 

for all a* < a ^ 1, the linear operator has a spectral gap of size More precisely, the 
spectrum of splits as &(Jfa) = {0} U 6 ^hh 

6 C {A G C ; Re A ^ — z/*} VaG(a*,l] (3.3) 

where we recall that Pq, denotes the spectral projection associated to the zero eigenvalue of 
tend denotes the part of on X = (1 — ¥a)X. 

To prove this we will use the following result asserting that if an operator has a spectral 
gap, and another operator is close to it in a certain sense, then it must also have a spectral 
gap of a comparable size: 

Lemma 3.4. Let X bea Banach space and let {Lq, S>{Lq)) he the generator of a Co-semigroup. 
For any e G (0,1) let {L^, ^(L^)) be a given closed unbounded operator with &{Lq) c S>{Lf) 
and 

lirn IKLg — Lo)-R(A,Lo)||^(x) = 0 VA G C with ReX > s{Lo). (3.4) 

Then, 

limsups(L£) ^ s{Lo) 

£—^0 

with 

lini ||ii(A, Lg) — ii(A, Lo)ll^(x) = VA G C with ReX > s{Lo). 

Proof Let A G C be given with Re A > s{Lo) and let sq > 0 be such that 

||(Lo ~ -^o)ll^(x) Ve G (0,eo). 

Setting 

= Lo)R{X, Lo) =:I-Z, 

one gets that is invertible for any e G (0, eo) with ^e- Since moreover 

/e = {X-L,)R{X,Lo) 
one gets that A — is invertible for any e G (0, eo) with 

R{X,L,) = R{X,Lo)/r^ VeG(0,eo). 

Therefore, Re A ^ s{Lf) for any e G (0,eo) which proves the first part of the result. For 
the second part, one has simply, for a given A G C with Re A > s{Lo), 

||i?(A,L,) -R(A,Lo)||^(^) = ||R(A,Lo) (A"' -^)IU(x) 

and it suffices to prove that 

j™IUr‘-z|Ux) = »' 

Since = Y.n=Q we get 

CXD 

~ ^ W^e ll^(X) 

n=l 

and clearly, since lim^^o ll^ell^(x) = 0 we get the conclusion. □ 
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Remark 3.5. We notice that the above result can also be seen as a way of stating general 
abstract results for relatively bounded operators (see [13, Theorem 3.17, p. 214]). 

We are now ready to prove Theorem 3.3: 


Proof of Theorem 3.3. We will apply Lemma 3.4 to the restriction of to X for a close 
to 1. (We recall the reader that the spaces X and y were defined in section 2.4.) Notice 
that, since 



du = 0 


V/ G 


a G ( 0 , 1 ] 


one can define the restriction .ifo,: d X ^ X hy Ci X = y and 

yfaf = yCaf for any / G T’ for any a G (0,1]. According to Theorem 2.6, < 0. 

Estimate (3.4) in Lemma 3.4 for and is exactly Proposition 3.1 since 


= \\yf^ih) - .^i{h)\\x v/i G 55. 


Since R{X,yfi): X ^ y, the hypotheses of Lemma 3.4 are satisfied and therefore s{.^a) ^ 
V* < V = s{^i) for any a close enough to 1. Now, since X = Range(I — Pq,) for any 
a G (oi, 1], one has = ■^a|(i-p„)A' for any a G (ai, 1]. This finishes the proof. □ 


3.2. Decay of the associated semigroup. Now, one should translate the above spectral 
gap of the operator into a decay of the associated semigroup. To do so, we use a stable 
splitting of the generator into a dissipative part and a regularising part. This strategy 
is inspired in the recent results [12, 18], but we give a proof adapted to our situation, 
exploiting a well known stability property of the essential spectrum under weakly compact 
perturbations. Our splitting is in the spirit of the one described in Section 2.4. Namely, set 

r«/ = QM, Fa) + Qa{F^, /), fey 

so that yfa = 7^ + L. The positive part of this operator is 

T+f = Qtif,F^) + Qf{F^,f) 

and Tq is written as 

Taf{v) = Tf'f{v)-Fa{v)[ f{w)\v-w\dw - aa{v)f{v), vdR^,fey, 

7r3 


where 

0'a(v) = / Fa{w)\v — w\dw y aa{l + \v\) Vu G 

Inspired by the splitting in Section 2.4, let us pick R> D large enough so that (2.13) holds 
true and define, for any a, 

Faf{v) = {XB-J){v) + T+{xb-J){v) - (s:('c) + (Ta{v)) f{v) 

for / G 3^ and G Let us first see that Ba thus defined is dissipative: 


(3.5) 
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Lemma 3.6. Let R > 0 and /3 > 0 be given as in (2.13). For any 0 < /3* < there exists 
(/3*) G (ao, 1) such that 


signf{v)Bafiv)m^{v)dv^-/34f\\y V/G 3^, VaG(a^l). 


(3.6) 


Proof. The proof is a direct consequence of (2.13) together with the fact that converges 
strongly to More precisely, let us fix / G T’ and compute first \\Ta fn — JrWx- One 
checks easily that 

\\r+fR-T+fR\\x ^ \\QtiFa-MjR) + QfifR,Fo,-M)\\x 

+ WQUMJr) - Qt{MjR)\\x + \\QtifR,M) - Qt{fR,M)\\x 

^ C \\Fa - M\\y \\fR\\y + 2p(l - a)||-Ad||^Li(^-i) ||/it||y 

where we used both Proposition 2.3 and 2.4. Consequently, there exists a nonnegative 
function (ii(a) with lima^.i (5 i(q:) = 0 such that 

\\Ta!^fR-T^^fR\\x ^Si{a)\\f\\y VaG(ao,l) Vi? > 0. (3.7) 

Set now 

= [ signf{v)Baf{v)m~^{v)dv VaG(ao,l]- 

iK3 

Using the fact that 

Baf{v) -Bif{v) = T^fR{v) -T^^fRiv) - (cFaiv) -cri{v))f{v) 

we get readily that 

^a(/) ^ ^l(/) + WT^fR - T+JrWx 

-[ (cTaiv) - ai{v)) \f{v)\m~^{v)dv 

iR3 

^ + 5i(a) WfWy + [ Wa{v) - o-i(i;)| \ f{v)\m~^{v) di; 

iK3 

where we used (3.7). Finally, since |r; — re] ^ (v) (w) Vu, w G M^, we have 

\aa{v) - Ui(i;)| ^ / \v-w\ \Fa{w) - Al(r(;)| dw ^ (r’)llF’a - Al||^im3) ^ {v)\\Fa - M\\y 

iR3 n J 

and we deduce from Theorem 2.5 that 

[ |cr«(i;) - CTi('(;)| |/(r;)|m“^(i;)di; ^ 7?i(a) f {v)\f{v)\m~^{v)dv = rii{a)\\f\\y 
iR3 iR3 

with limQ^iT/i(a) = 0. To summarize, there exists a function 6{-) with limc^i <5(a) = 0 
such that 

^a(/)^^i(/) + i^(a)||/b V/g3^ 

which, from (2.13), becomes 

W«)-/3) WfWy V/g3^. 
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This gives the first part of result since lima^-i 5(a) = 0. □ 

We set now Aa = -^a — ot G (al, 1), or in other words 

•^o^f{v)=T^{XBRf){v)+l^'^{XBRf){v)-Fa{v) [ f{w)\v -w\dw 

Jm.3 

for u G and f G A. Then we have then the following: 

Proposition 3.7. For any a* G (0, o) and for any a G {of, 1), one has 

i) Aa G ^{X); 

ii) Ba : S>{Ba) C X ^ A with domain ^{Bq) = y is the generator of a Co-semigroup 
{Uait))t^Q of A with 

Wa{t)f\\x ^ex.p{-f3A)\\f\\x Vt ^ 0, V/G Af. (3.8) 

Proof It is clear from Proposition 2.3 that Aa is a bounded operator in A since 

II^q/IIa: ^ CWFaWyWxBnfWy + l|-^a||yll/llLl(R3) 

^CR\\Fa\\y\\f\\x +\\Fa\\y\\f\\x '^fGA 

for some positive constant Cr depending on R. 

Since SAa (with domain y) is the generator of a Co-semigroup in A according to The¬ 
orem B.3 and Aa is bounded, it follows from the classical bounded perturbation theorem 
that Ba (with domain T’) is also the generator of a Co-semigroup {Ua{t))^^Q in A. Since 
Ba + (d-k is dissipative according to (3.6), (3.8) holds according to the Lumer-Phillips theo¬ 
rem. □ 

Actually, it is easy to check that Aa has better regularising properties: 

Lemma 3.8. For any a G (ao, 1), Aa G idS{A,y). Moreover, there exists some Maxwellian 
distribution M such that, for any a G (ao, 1), 

Aa G ^{A, n) where R = 

In particular, Aa is a weakly compact operator in A. 

Proof Recall that there exist two Maxwellian distributions M and M. (independent of a) 
such that 

M(v) ^ Fa{v) ^ M.{v) \/v G Va G (ao, 1). 

In particular, there exists some Maxwellian distribution M such that 

sup \\Fa\\a(M-'^)= sup ([ M~^{v){vf \f{v)\‘^ dv] =Cm<oo. 

a6(oo,l) «S(ao,l) / 
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Then, because |i; — l(;| ^ (v) (w) for any v,w ^ one has first that, for all / E Af, 


/ Fa{v) / f{w)\v — w\dw 

i]R3 J]R3 






\Fa{v)f{v)‘^M^{v)dv f I \f{w)\{w) dw 


^ WFn, 




Moreover, according to [1, Proposition 11], there exists C > 0 such that 

||Q^(5,^)|L2 (m-i) ^ and 

Using this with g = fxBu h = FaWe get that there exists C > 0 (independent of a) 
such that 

\\-^aifXBR)\\j^2(^M-^) ^ ^ (ll /XSflllLi(M-i/2) + ll/XSflllLi(M-i/2)^ 

In particular, there exists C = Cr > 0 such that 


In the same way, 

^ CrII/iI;^ v/ea”. 

This proves that Aa E Due to the Dunford-Pettis theorem the embedding % 

X is weakly compact, which proves the second part of the Lemma. □ 


We this in hands, one has the following result about the decay of the semigroup 
{Sa{t))i-^Q in X generated by Ffa- Remember that u is the spectral gap of the elastic 
linearised operator jSfi. 


Theorem 3.9. Take 0 < Z 2 * < Z 2 and 0 < /3* < a (where (3 > 0 is such that (2.13) holds). 
With the notations of Theorem 3.3 and Lemma 3.6, let ai = max(a*, a'*'). Then, for any 
a E (ai, 1) and for any g E (i^*, z^), there exists C = C{g, a) > 0 such that 

\\SUt){I-Fa)\\^ix)^Ce-^' Vt^O 

where Pq, is the projection operator over Span(FQ,) in X. In other words, for any ho G X the 
solution h = h{t, v) (in the sense of semigroups) of the equation 

dth = 


satisfies 


||/i(t) - cG^IIa’^ C'll/iolle fort^O, 


where c := /r 3 ^o- 
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Proof. Denote by {Ua{t))t^o the semigroup in X generated by Ba- Since = Aa + Ba, it 
is well known from Duhamel’s formula that 

Sa{t) = Uait) + f Sa{t - s)Aalia{s)ds. 

Jo 

Since Aa is weakly compact in X, one gets that, for any t ^ s ^ 0, the integrand 
Sa{t — s)AaUa{s) is a Weakly compact operator in X. Using then the “strong compact¬ 
ness property” (see [9, Theorem C.7] for general reference and [19] for the extension to 
weakly compact operators in L^-spaces) we get that 

Sa{t) — Uait) is a weakly compact operator in X for all t ^ 0. 

We recall that, by definition, the Schechter essential spectrum is stable under compact 
perturbations. However, it can be shown that in L^-spaces it is actually stable under 
weakly compact perturbations, see [14, Theorem 3.2 & Remark 3.3]. Due to this property 
we have 

®ess (5a(t)) = ®ess (Uait)) Vt ^ 0. 

In particular, the two Co-semigroups share the same essential type, i.e. uiessiSa) = (^essiUa)- 
Since coessiUa) ^ wo(Z^a) we get 


^688(^50) ^ / 3 * < 0 . 


Since 


9^0 (5q,) = max(u;es8(5Q,),s(.ifa)) 
with s(Jfa) = 0 we obtain that 


^0(‘5q) - 0 > ^@88(^50.). 

General theory of Co-semigroups [9, Theorem V.3.1, page 329] ensures that, for any w > 
^ess(^a)) one has 

6(jSf„) n {A G C ; Re A ^ Cj} = {Ai,..., A4 

with Aj eigenvalue of Jfa with finite algebraic multiplicities ki and Re Ai ^ ... ^ Re A^ 
and there is C^; > 0 such that 

||5Q,(t)(/- n^)|| ^ C^ exp(a;t) Vt ^ 0 

where is the spectral projection associated to the set {Ai,..., A^}. In particular, choos¬ 
ing 1 = 2, since Ai = 0 while Re A 2 is the spectral gap of (see Theorem 3.3), choosing 
then 0 < p < 

6(if„) n {A G C; Re A ^ -/r} = {Ai} = {0} 

we get the result since = Pq, is the spectral projection on the simple eigenvalue 0. □ 
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4. Local stability of the steady state 

Using the result of the previous section, we show here that when a is close to 1 and 
the initial condition is close to the steady state, solutions to Eq. (2.8) converge to it 
exponentially fast. Choosing a close enough to 1 and the initial condition close enough to 
equilibrium, the exponential speed of convergence can be as close to v as we want. In all 
this section, we shall assume that the initial datum /o is such that there exist b > 0,s e 
( 0 , 1 ) satisfying 


/ /o(r^) exp( 6 |r;|) df < oo. (4.1) 

Jms 

Then, according to Theorem 2.1, there exists 0 < a < 6 (independent of a) such that the 
solution f(t) of (2.8) satisfies 

f(t) G T” Vf ^ 0 

where we recall that 

m(i;) = exp (—a|i;l), (o > 0 ), and = L^(U^,m~^(v)dv). 

Actually, assuming a bit more on the initial datum, one can prove the following where we 
recall that y = L\{m~^)\ 

Lemma 4.1. Let f{t, v) be the solution to (2.8) associated to a nonnegative initial condition 
fo satisfying (1.4) and (4.1). Then, for any e > 0 there exists C > 0 depending only on e 
and the moment (4.1) of fo such that 

\\fmy^C\\f{t)\\]f^ Vf^O. (4.2) 

Proof By Holder’s inequality, taking a dual pair p, q, i.e. l/p + l/q = 1 


\\fit)\\y= f{t,v){v)m (i;)di; 




\ i/p 

f{t, v)m~^{v) dt! j 


\ !/<? 

f{t, v){vYm~^{v) dv j 




The initial datum fo has an exponential tail of order 1 and Theorem 2.1 ensures that this 
property propagates with time; hence for any p ^ 1 there exists C > 0 such that 


sup 


\ 1/9 

f{t^v){vYm~^{v)dv\ ^ C < oo, 


where the constant C > 0 depends only on q and fo (and not on a). Hence, choosing p 
such that 1 /p > 1 — e proves the lemma. □ 


We deduce from the above lemma the following stability result: 
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Theorem 4.2. Take 0 < i/* < i/. There exist 0 < a* < 1 and e > 0 such that, for all 
a with a* ^ a < 1 and all nonnegative /o € 3^ satisfying (1.4) and (4.1) and such that 
Wfo-Fj X ^ die solution f{t) of (2.8) with initial data /o satisfies 

\\f{t) - Fa\\x ^Cexp{-u^t)\\fo-Fa\\x (t ^ 0 ) (4.3) 

for some constant C which depends only on a and the moment (4.1) of the initial condition 
fo. 

Proof Let a G (ao, 1] and let f{t, v) be the solution to (2.8) associated to the initial datum 
ft). Setting h{t,v) := f{t,v) — Fa{v), we may rewrite (2.8) as 

dfh = + Qa{h, h), 

and through Duhamel’s formula, denoting by (5o(t))^^o the semigroup generated by Jtfa, 

h{t) = Sa{t)[h{0)] + f Sait - s)[Qa{h{s),h{s))] ds. 

Jo 

Take any p, such that v,, < p < v and consider oq < «* < 1 such that for every a* ^ 
a ^ 1 , the operator has a spectral gap of size p (the existence of such a* is warranted 
by Theorem 3.3). Then, for this range of a, the semigroup ((I — Pa) 5 a(f))i 3 :o decays 
exponentially with speed p. Recalling that Range(I — Pq,) = X and that h{t) € X for any 
t ^ 0 we get that there exists C > 0 which depends only on a such that, 

ll^(i)II.V ^ ll‘5«W[^(0)]|U + / ||5a(t- s)[Qa(/l(s),/l(s))]||_^ ds 

Jo 

^ C||/i(0)||_^ exp(-/it) + C [ \\Qaih{s),his))\\^ exp{-p{t - s)) ds 

Jo 

^ C\\hi 0 )\\^ exp{-pt) + C l|/i('S)||J exp(-^(t - s)) ds 

^ C'||/i( 0 )||_^ exp(-^t) + C's ^ ll^('S)||^^ exp(-/i(t - s)) ds, 

where we have used Lemma 4.1 with e = 1/4. 

From this point, a Gronwall argument is enough to show that for || /i(0) ||^ small enough, 
h{t) converges exponentially fast to 0, with a speed as close to p as we want. Let us develop 
this argument more precisely. Take (5 > 0 and fo such that ||/i(0)||^ ^ 6/2, and consider t 
in a time interval [ 0 ,T, 5 ] where ||/i(t)||^ ^ 6 . Then, 

\\Ht)\\x ^ C\\hi 0 )\\^ expi-pt) + 6 ^I’^C 2 f \\his)\\^ exp{-p{t - s)) ds 

Jo 

or, rewriting this for the quantity 7 (t) := exp(^t)||/i(t)||_^, 

7(t) ^ C7(0) + fl{s)ds (fG[ 0 ,r 5 ]). 

Jo 
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Then, by Gronwall’s Lemma, 

7(t) ^ 6*7(0) exp((5^/^C'2t) (t G [OjT^]), 


or equivalently. 


\\h{t)\\^ ^ C\\h{Q)\\^ exp(((5^/2C'2-/r)t) {t G [0,r5]). 

Now, take 6 small enough so that 5^/'^C2 — Ai < and e < (5/2 small so that 

Ce exp(— ^ 5. 


(4.4) 


Then, for ||/i(0)||^ < e, one can actually take = +oo and (4.4) finishes the proof. □ 


5. Global stability 


5.1. Evolution of the relative entropy. We consider the evolution of the relative entropy 
of a solution /(t, v) to (1.1) with respect to the equilibrium Maxwellian . Notice that M. 
is not the steady solution associated to (1.1) and f{t, v) is not expected to converge towards 
M. However, we shall take advantage of the entropy-entropy production estimate satisfied 
by L to estimate first the distance from f(t) to A4 which, combined with Theorem 2.5, 
yields a control of the distance between f(t) and Fa- Using the entropy-entropy production 
estimate. Theorem 1.3, obtained in [5], one has the following crucial estimate on the 
evolution of the relative entropy along solutions to (1.1) 

Proposition 5.1. For any a G (oq) 1] CLnd any nonnegative initial datum /o with unit mass 
and H{fo\M) < oo, the solution f{t) = f{t,v) to (1.1) satisfies: 


H(f{t)\M) < exp(-At)iT(/o|Al) + K(1 - a) Vf ^ 0 ; a G (oq, 1] 


for some positive constant K > 0 independent of a where X > 0 is the constant appearing in 
Theorem 1.3 

Proof Let a G (ao, 1] be fixed. Given a solution f{t) = f{t, v) to (1.1), set then 



Gomputing the time derivative of H (t) we get 



We recall [11] that 
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where the entropy production functional &H,a{g) is defined, for any nonnegative g, by 


^H,a{9) = -^ [ 

oTT jRa 


\v - w\g{v)g{w) 

' g{v')g{w') 


log5KMKl-lldud.da;>0 


g{v)g{w) g{v)g{w) 

where the post-collisional velocities {v',w') = are defined in (2.4). Moreover, 

using the definition of M. and the fact that Qa conserves mass, one has 

1 


Qo,{fJ)\ogMdv = - 


= +- 


20 Jm 

1 — a" 
160 


Qa{fJ)\v\ du 


■ / f{t,v)f{t,w)\v — w\^ dvdw 


(5.2) 


where we used (2.2)-(2.3) noticing that Aa [| • P] {v,'w) = — ^ |f — tcp. Putting (5.1) 
and (5.2) together we obtain 


= -D(f(t}) - ^H,a(f(i)) + ^o 2 [ f{t,v)f{t,w)\v-w\dvdw 

dt 2a^ JrSxrs 

1 ^ f 

-—^ / f{t,v)f{t,w)\v - wf dvdw 

16c) JrSxrs 

1-0 f 

-^ / f{t,v)f{t,w)\v — w\dvdw. 

2(^0 ^R3xR3 

Using the uniform control of moments of / one sees that there exists some positive con¬ 
stant C > 0 independent of a such that 

^ -XH{t) + C{l-a) Vt ^ 0 , Vo G (oq, !]• 

Integrating this inequality, we obtain 

H{f{t)\M) ^ exp{-Xt)H{fo\M) + (1 - exp(-At)) 

which gives the result. □ 

5.2. Global stability result. Using the Csiszar-Kullback inequality we deduce from Propo¬ 
sition 5.1 the following 

Theorem 5.2. There exist some T = T{ao) and some function I : (aojl] — with 
lima^i i{a) = 0 such that 

\\f{t) - FaWx ^ -((a) Vt ^ r(ao); « S (oq, !]■ 

We remark that both T and I in the above theorem can be given explicitly. 
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Proof. Using both the Csiszar-Kullback inequality and Holder’s inequality we get 

Due to Theorem 2.1 (and recalling the definition of m by (4)), assuming a <rl2 gives 

sup ||/(t) - MWilUn-^) < OO. 

t^O 

Using now Proposition 5.1 we get the existence of two positive constants Ci, C 2 > 0 
independent of a E {ao, 1] such that 

\\f{t) -M\\x ^ C'iexp(-Af) + C' 2 (l - a) Vt ^ 0; a E (ao, !]■ 

This, combined with Theorem 2.5, yields 

\\f{t) - Fa\\x ^ Ciex.p{-Xt)+ r]i{a) Vt ^ 0 ; a E (oq, 1] (5.3) 

where rjif) is a given explicit function with lima^i r/i(a) = 0. We get readily the conclu¬ 
sion. □ 

The previous results essentially contain the proof of our main result; 

Proof of Theorem 1.2. With the notation of Theorem 4.2, one can pick ai E (ao, 1) so that 
£(a) ^ e for any a E (ai, 1) so that 

\\f{t) - FaWx ^ e Vt ^ r(ai) : Va E (ai, 1] 

and Theorem 4.2 yields the global stability result. □ 
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Appendix A. Some results in perturbation theory of linear operators 

We gather here some results that are needed in Section 3.1 in order to study the spectral 
properties of Ff a for a close to 1. We begin by defining the gap between two closed linear 
operators, following [13, IV.2.4, p. 201]: 

Definition A.I. Let X, Y be Banach spaces and S, T closed linear operators from X to Y. 
Let G{S), G{T) be their graphs, which are closed linear subspaces of X xY. We set 

5{S,T) ■.= 5{G{S),G{T)) ■.= sup dist(n, G(r)), 

u&G{S) 

and we define the gap between S and T as its symmetrisation: 

6{S, T) := max{,5(5, T),5{T, 5)}. 

We include here Theorem 2.14 from page 203 of [13] for the convenience of the reader; 
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Theorem A.2 ([13, Thm. 2.14, p. 203]). Let X, Y he Banach spaces and A, T be closed 
operators between X and Y such that A is T-bounded with relative bound less than one; that 
is, 

||ALt||y ^ u G &{T). 

for some o ^ 0, 0 ^ 6 < 1. Then T + Ais a closed operator and 


5{S,T)^ 


yj a? + 1)^ 
l-b 


Finally, the following theorem is the main perturbation result we use in order to deduce 
the properties of the spectrum of Jtfa' 


Theorem A.3 ([13, Thm. 3.16, p. 212]). Let T he a closed linear operator on a Banach 
space X and assume its spectrum &{T) is separated into two parts by a closed curve F in C. 
Let X = X'rp 0 X'f he the associated decomposition of X. Then there exists <5 > 0, depending 
on T and T, such that any operator S on X with 

6(8, T) < 6 

satisfies the following properties: 

(1) The spectrum &(S) is also separated into two parts by the curve F. 

(2) In the associated decomposition X = XgtBX'f, the spaces X( and X'^ are respectively 
isomorphic to Xf and X'f. 

(3) The decomposition X = X^ © X'f is continuous in S in the sense that the projection 
Ps of X onto X'g along X'f tends to Pt in norm as 5(8, T) —)■ 0. 


Appendix B. Proof that generates an evolution semigroup 


Let 

be fixed and let 


m(v) = exp(—a|v|), a > 0 


X = L^(m ^(?;)di;), y = L^({v)m ^(v)dv). 

We wish here to investigate the compactness properties of the ’gain’ part of Afa for a G 
(ao, 1], with the final aim of showing that yfa generates a semigroup for all 0 < a ^ 1. 
Recall that 

^ah = Qa(h,F^) + Q^(F^,h)+L(h), hey. 

and also that 

L(/i) =/C/i-SO 

where 

K,h(v) = Q^(h,M)(v) = / k(v,w)h(w)dw and S(?;) = / Al(t(;)|v — t(;| dtc 

with 


k(v, w) = Cf\v — tel ^ exp I —/3o I |v — rc| + 


I |2 I |2 

V — u\ — \w — u\ 


\V — W\ 


(B.l) 
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with /3o = ^ and Co > 0 a positive constant (depending only on ©o). Notice moreover 
that S(r;) = k{w, v) drc. In the same way, one can write 

Qa{h,^a) = K^ah-aa{-)h and Qa{Fa,h) = IClh - IC^h 

with 

ICih=Qt{h,Fa), JClih) = QtiF^,h) 

while 

(^aiv) = / FQ,(t(;)|r; — rc| drc and }C^{h){v) = Fa{v) / h{w)\v — w\dw. 

With this notation 

^ah = lCh + K\h + lClh-{^ + aa)h-lClh 'ih^y. 

As for L, the two operators i = 1,2 are integral operators with explicit kernels. 
Namely, 

Lemma B.l. For any h Gy, one has 

}C^h{v) = [ K^{v,w)h{w) dw 

jrs 

where 

kI^{v,w) = I I [ Faiv + V 2 + ^ dV '2 (B.2) 

P — w\ Jv2-(w-v)=o V a + 1 J 

for some positive constant Ca > 0. 

Proof The proof follows standard computations performed for instance in [2] where 
was replaced by a given Maxwellian. In particular, (B.2) is derived in [2, p. 524]. 

Recalling that there exist two Maxwellian distributions M_ and M. (independent of a) 
such that 

M_{v) ^ Fa{v) ^ M.{v) Vf G Va G (0,1). 

In particular, this proves that, for any h y 0, 

}Cih^lclh=Qf{h,M) 

and 

)Clh^lclh = Qf(M,h). 

Again, JC^ h is an integral kernel with explicit kernel, namely 

{v, w)h{w) dtc 


with 

k\^{v, w) = Ca\v — exp < —/3i 


(1 + p,a)\v -w\ + 


\V — Ul\ 


\W — Ul 


\v — w\ 



JCihiv)= [ kI 

JR3 


(B.3) 


□ 0 
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where > 0 is a positive constant depending only on a and M while 

1 


1 — a 

/ia = 2^— ^0, /3i = 


'1 + a " "" 801 

01, ui being the kinetic energy and momentum of M. (see [2]). By a simple domination 
argument (namely, Dunford-Pettis criterion), if are weakly compact in X then so will 
be Ki^, i = 1,2. 


Proposition B.2. Let a G (ao, 1) be fixed. Then, 

1C y^x, ici : y^x 

are positive, bounded, weakly compact operators. Moreover, G B{X) while 

JCl : ^ 

is a bounded and weakly compact operator. 

Proof The fact that 1C, Kf are bounded operator from y to X comes from Prop. 2.3. We 
divide the proof of the compactness properties into several steps. 


First step: weak compactness of Kf. We already notice that it is enough to prove that 
ic\ : y ^ X is weakly compact. Let dv{v) = m~^ du and let B = By be the unit ball 
of T’. Since X = L^(R^, dv), according to Dunford-Pettis Theorem, this amounts to prove 
that 


sup f }Cl^h{v) di>{v) —)■ 0 
h&B Ja 


as h>{A) -A 0 


and 


sup / 

J\v—ui\>r 


lC^h{v) dvfv) —>• 0 as r ^ oo. 


(B.4) 


(B.5) 


Using the representation of JC^ as an integral operator, it is easy to check that (B.4) and 
(B.5) will follow if one is able prove that 


and 


sup [ K]^{v,w)m ^(r;)dr'—^0 as v{A) ^ d 

ioeR3 \^) J A 

m{w) f 


-ttI 


sup . 

uigR3 \^/ J\v—ui\>r 


K^{v,w)m ^(i;)dt;—>■ 0 as r —>■ oo. 


(B.6) 


(B.7) 


Let us prove (B.6). Let A c be a given Borel subset and let tc G be fixed. Set 

Bw = {r” G , \v — w\ < 1}. Since k\{v, w) ^ Ca\v — one has 


/ K\{v,w)m ^ (v) df ^ Cq, / \v — w\ ^ du{v) 

Ja Ja 


= Cn 


' Ar\B,, 


V — w\ ^ dv{v) + 


'Ar\B<r„ 


t; — u;| ^ du{v) . 
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Clearly 


/ \v — w\ ^ 011 ^( 7 ;) ^ i^(^) 


while, for any p > 1, 1/g + 1/p = 1, one has 


' AnBu 


V — w\ ^ ^ 


'AnBu 


di'{v] 


1/9 


'AnB^ 


\v — wl ^ exp(a|t'|) 


i/p 




exp 




— w\ ^ dv 


i/p 


where we used that, exp(a|?;|) ^ exp(a|r(;|) exp(a|r' —wl) for anyUjW. Choosing nowp > 3, 
one sees that 

\ i/p 


V — w 


-p 


dr; ) < 


00 


' B^ J 

and is independent of w. Thus, there exists C = C{a, a,p) such that 


kI 


{v,w)m ^(f) dv ^ C ^exp Vw E 


Since p > 1, this proves that (B.6) holds true. Let us now prove (B.7). One first notice 
that 


K‘^{v,w)m ^(t;)d'(;—)• 0 as r — 00 . 


miw) f 
sup / 

|io—ttl|s;r/2 \^/ J\v—ui\>r 

Indeed, one notices that 

kI^{v,w) ^Ca\v - w|“^ exp (2/3i(l + Pa) [|w - -\v - Mlp]) • 

Therefore, if |w — rii| ^ r/2 and |ri — rti| > r, one gets 

k\{v,w) ^ —^exp (-|A (1 + ^J-a)\v - Ui\^) 

and (B.8) follows easily since 

/ exp (—|/3i(l + Pa)!?; — Mi|^) di^(r') < 00 . 

Jm? 

Now, to prove (B.7), it is enough to show that 

m{w) f 


(B.8) 


sup 

\w—u\\>r/2 \^) 


K‘^{v,w)m ^(?;)d?;—)• 0 as r —>■ 00 . 


(B.9) 


' \v—ui\>r 

Arguing as in [4, Proposition A.l] (with s = 1), there exists K = > 0 such that 


K^{v,w)m '^{v)dv^Kam (w). 


Therefore, for any r > 0, 

m{w) 


i"^) J\v—ui\>r 


K^{v,w)m '^{v) dv ^ Ka{w) 


-1 
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and (B.9) follows since sup|^_„j |>^/2 (^) ^ 0 as r —oo. This achieves to prove that 

• 3^ —^ -T is weakly compact. 


Second step: weak compactness of 1C. Notice that K, and /C^ are two integral operators 
whose kernels, given respectively by (B.l) and (B.3), are very similar. The above compu¬ 
tations can then be reproduced mutatis mutandis to get the weak-compactness of K,. 


Third step: boundedness of IC^. According to [1, Theorem 12], and since m{v) = exp(—alul^) 
with s = 1, one has 

K = Qt(Mr) : 

is bounded for any 0 < a < 1. Consequently, simple domination argument asserts that 

ICl E B{X). 


Final step: weak compactness ofKf^. Recall that 

K?^h{v) = ¥a{v) / h{w)\v — w\dw. 

Therefore, 

\IClh{v)\ < (v)Fa(u) [ |/i(rc)|(u;)dr(;. 

In particular, there exists C > 0 such that 

< {v)Fa{v) [ \h{w)\m~^(w)dw. 

Since 

/ FQ,(u)(ti)m“^(r;) dt; < oo 

this proves that JC^ : X ^ X is bounded and dominated by a one-rank operator. In 
particular, it is weakly compact. □ 

As a general consequence, one has the following 


Theorem B.3. For any a E (ao) 1)^ the unbounded operator 5^^ is the generator of a Co¬ 
semigroup {Sa{t))^^Qin X. 

Proof One applies the recent version of Desch theorem for positive semigroups in L^- 
spaces, see for instance [20]. For any a E (ao, l)j define the multiplication operator: 

Aa h{v) = -{cFa{v) + T,{v))h{v), h E ^{Aa) = y 

then, Aa is the generator of a positive Co-semigroup {Ua{t))t^o in X. Since K, and Kf are 
weakly compact, one has from [20] that Ba = Aa + K, + Kf is the generator of a positive 
Co-semigroup {Va{t))t^o- Finally, since and are bounded operators, one gets that 

.^a = Ba +ICl- ICl 

is the generator of a Co-semigroup in X. □ 


STABILITY OF THE STEADY STATE 


30 


References 

[1] Alonso, R. J., Carneiro, E. & Gamba, I. M. Convolution inequalities for the Boltzmann collision 
operator. Com. Math. Phys., 298 (2010), 293-322. 

[2] Arlotti, L. & LODS, B., Integral representation of the linear Boltzmann operator for granular gas 
dynamics with applications. J. Statist. Phys. 129 (2007), 517-536. 

[3] Bisi, M., Carrillo, J. A. & Lods, B., Equilibrium solution to the inelastic Boltzmann equation driven 
by a particle bath, J. Stat. Phys. 133 (2008), 841-870. 

[4] Bisi, M., Canizo, j. a. & Lods, B., Uniqueness in the weakly inelastic regime of the equilibrium state 
to the Boltzmann equation driven by a particle bath, SIAM Journal of Mathematical Analysis, 43 (2011), 
2640-2674. 

[5] Bisi, M., Canizo, J. A. & Lods, B., Entropy dissipation estimates for the linear Boltzmann operator, to 
appear in Journal of Functional Analysis, 2015. 

[6] Brilliantov, N. V. & Poschel, T., Kinetic theory of granular gases, Oxford University Press, 2004. 

[7] Carrillo, J. A. & Toscani, G., Contractive probability metrics and asymptotic behavior of dissipative 
kinetic equations, Riv. Mat. Univ. Parma 1 (2007) 75-198. 

[8] Edmunds, D. E. & Evans, W. D., Spectral theory and differential operators, Oxford University Press, 
New York, 1987. 

[9] Engel, K. J. & Nagel, R., One parameters semigroups for linear evolution equations. Springer, 
1999. 

[10] Ernst, M. H. & Brito, R., Scaling solutions of inelastic Boltzmann equations with over-populated high 
energy tails. J. Stat. Phys. 109 (2002), 407-432. 

[11] Gamba, L, Panferov, V. & Villani, C. On the Boltzmann equation for diffusively excited granular 
media. Comm. Math. Phys. 246 (2004), 503-541. 

[12] Gualdani, M.P., Mischler, S. & Mouhot, M. Factorization for non-symmetric operators and expo¬ 
nential Jf-theorem, preprint 2011, http://hal.archives-ouvertes.fr/hal-00495786/en/. 

[13] Kato, T., Perturhation theory for linear operators. Second edition. Classics in Mathematics, Springer 
Verlag, 1980. 

[14] Latrach, K., & Jeribi, a.. Some Results on Fredholm Operators, Essential Spectra, and Application, J. 
Math. Anal. Appl. 255 (1998), 461-485. 

[15] Mischler, S. & Mouhot, C., Cooling process for inelastic Boltzmann equations for hard spheres. 
Part II: Self-similar solution and tail behavior, J. Statist. Phys. 124 (2006), 703-746. 

[16] Mischler, S. & Mouhot, C., Stability, convergence to self-similarity and elastic limit for the Boltzmann 
equation for inelastic hard spheres. Comm. Math. Phys. 288 (2009), 431-502. 

[17] Mischler, S. & Mouhot, C., Stability, convergence to the steady state and elastic limit for the Boltz¬ 
mann equation for diffusively excited granular media. Discrete Contin. Dyn. Syst. 24 (2009) 159-185. 

[18] Mischler, S. & Scher, J., Spectral analysis of semigroups and growth-fragmentation equations. Ann. 

I. H. Poincare -AN (2015), to appear, http://dx.doi.Org/10.1016/j.anihpc.2015.01.007 . 

[19] Mokhtar-Kharroubi, M. On the convex compactness property for the strong operator topology and 
related topics. Math. Methods Appl. Sci. TJ (2004) 687-701. 

[20] Mokhtar-Kharroubi, M., New generation theorems in transport theory. Afr. Mat. 22 (2011) 153-176. 

[21] VAN Noije, T. & Ernst, M., Velocity distributions in homogeneously cooling and heated granular fluids. 
Gran. Matt. 1 (1998), 57-64. 

Jose A. Canizo, Departamento de Matematica Aplicada, Universidad de Granada, Av. Fuentenueva 

S/N, 18071 Granada, Spain. 

E-mail address: canizo@ugr.es 

Bertrand Lods, Dipartimento di Statistica e Matematica Applicata & Collegio Carlo Alberto, 

UniversitA degli Studi di Torino, Corso Unione Sovietica, 218/bis, 10134 Torino, Italy. 

E-mail address: lods@econ.unito.it 



